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Abstract 

This paper deals with the problem of oscillating flows occurring in devices such as Stirling or thermoacoustic engines and 
refrigerators. Since the global governing equations cannot be solved, the authors propose to introduce a few simplifications; the 
most simplifying reduction is that the fluid is assumed to be incompressible. However, specific attention is paid to describing 
the flow characteristics that’s why the Lagrangian formalism which allows the individual study of each fluid particle is adopted. 
Thereby each particle contribution to global thermal effects can be evaluated and the gas temperature profiles along the exchanger 
can be computed. Various situations are presented including the case of a non-uniform temperature at the wall and a phase lag 
between pressure and temperature at the fluid entrance. The efficiency of the wall to fluid thermal exchange is analyzed. The 
authors show that this exchange depends upon two important parameters: the geometric ratio between the exchanger length and 
the particle oscillating displacement, and a thermal parameter “/3”, governing the temperature profiles and related to the Prandtl 
number, the operating frequency and the phase lag between the instantaneous heat flux and the wall to fluid temperature difference. 
© 2007 Elsevier Ltd and HR. All rights reserved. 
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Nomenclature 

c specific heat capacity at constant pressure 

(Jkg-'K- 1 ) 

Cf friction factor 

D gas displacement (m) 

d h hydraulic diameter (m) 

e thickness (m) 

E exchange effectiveness 

F\. F\ functions Eqs. (22) and (34) 

H convection heat coefficient (W m~ 2 KA 1 ) 

h,H local and radial average of enthalpy flux 

(W m~ 2 ) 

j complex imaginary (j 2 = — 1) 

/ length (m) 

Nu, Nu r i. Nil/ j Nusselt number, real part, imaginary part 
P pressure (Pa) 

Pr Prandtl number 

Q thermal energy (J) 

Re Reynolds usual number 

Spass internal cross-section (m 2 ) 

T time (s) 

T temperature (K) 

t* = tot non-dimensional time 

u, v velocities for x and y directions (m s _1 ) 

u* = (2 u/coD) non-dimensional fluid velocity 

fPom Womersley/Stokes number 

x axial coordinate (m) 

x* = (x/(lo/2)) non-dimensional axial coordinate 

y, z co-ordinates (m) 

y* = (y/(e/ 2)) non-dimensional co-ordinates 


Greek 

X parameter Eq. (11) 

Xj thermal conductivity i (W m 1 KA 1 ) 

H dynamic viscosity (kg m _1 s _1 ) 

p density (kg nA 3 ) 

w pulsation (rad s -1 ) 

f Lagrangian variable 

A thermal gradient (K nA 1 ) 

AT rad radial gradient T g — T i (K) 

APiong longitudinal gradient T\ — 7A \ (K) 

Ml fluid relative displacement 

Subscripts and superscripts 
8 gas 

in, in' entrance in x* = — 1 and in x* = 1 

/o extremity in x = l 0 

rad, long radial and longitudinal 
t tube 

0 initial, for x = 0, of reference 

out, out' Exiting in x* = — 1 and in x* = 1 
1,-1 to right or left extremity 

Operators 

p*, T* normalized magnitude with respect to reference 
value p 0 , T q 

3 imaginary part 

K real part 

x x cross-section average 

|a'| x magnitude 

(x) x temporal average on a semi-period 


1. Introduction 

Numerous varieties of thermal machines operate with 
an internal oscillating fluid flow: Stirling cooler [1], all 
type of Stirling engines including the kinematic engines 
and the free piston or free displacer engines: Martini 
model [2], Ringbom model [3] and hot air flow engine 
like Manson model. All versions of Pulse Tube Refriger¬ 
ators (PTR, OPTR, DIOPTR [4,5], etc.) and thermoacous¬ 
tic systems based on standing or progressive waves 
(Ceperley ring [6]) are also concerned. There are also 
a lot of hybrid systems like Gifford-Mc Mahon and 
Vuilleumier coolers [7] or fluidynes pumps and liquid pis¬ 
ton Stirling machines [8,9]. All these machines are known 
for many years but are not of a common use yet. How¬ 
ever, they may take a bigger place in a near future 
because of their abilities to operate with various energy 
sources brought by external ways (hence removing most 
technical problems occurring in internal combustions): bio¬ 
gas, solar energy, geothermic, heat collected downstream 


industrial processes, and thermoacoustic waste heat 
upgrading. 

In most cases, the mechanical conception of these 
machines is quite simple compared to gas turbines, internal 
combustion engines or classic refrigerating and gas liquefac¬ 
tion devices. Note that some of them (purely thermoacoustic 
devices) do not require any mobile part. In this paper the 
authors propose to focus their interest on a particular part 
of these systems: the heat exchangers where at least one fluid 
operates in oscillating flow mode. Moreover, such devices 
can improve many industrial processes; for example numer¬ 
ous works [10,11] show the interest of oscillating flows in 
micro-heat exchangers used as coolers for electronic 
components. 

The oscillating fluid flow causes thermal exchange phe¬ 
nomena which are not well estimated yet. Indeed, due to 
a lack of knowledge, most researchers and engineers use 
classic heat exchangers calculus, assuming fluids moving 
in a steady state flow regime which is not in connection 
with reality. As mentioned above, very few theoretical or 
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experimental studies concerning oscillating flows in heat ex¬ 
changers have been developed during the past decades. 
That’s the reason why Gedeon [12] and Swift [13] can be 
considered as precursors. More recently Wakeland and 
Keolian [14—16] propose some interesting results on ther¬ 
mal transfers between two exchangers in a standing wave 
acoustic field. Similarly, in a precedent paper [17,18], the au¬ 
thors point out two fundamental aspects for designing such 
devices: first the exchangers length in oscillating flow 
mode must be equal to the fluid peak to peak displacement 
and second, the temperature gradient at the wall may intro¬ 
duce a “regenerator effect”. In the latter consideration, 
a fluid particle may exchange with the wall a heat quantity 
exactly equal to the heat quantity it will reload at the same 
location when flowing in the opposite direction, hence the 
global thermal effect is null. In fact, the main parameters 
of interest for describing laminar and quasi incompressible 
flows are the fluid relative motion A L , which depends of 
the mechanic configuration and the Womersley number 
W om , which is linked to the operating frequency [17]. The 
Womersley number represents in fact the ratio of inertial 
forces to viscous forces in non-stationary regime, it is also 
an information on the relative thickness of dynamical 
viscous boundary layers compared to hydraulic radius of 
the canal. 


W = 

yr om 



(i) 


In addition complex forms of Nusselt numbers are used to 
characterize the heat transfers between walls and oscillating 
fluids. This allows the introduction of the phase lag between 
temperature variation and wall thermal exchange. 

The present paper is an extension of Ref. [17]: after 
a brief description of the most significant modeling steps 
and general results, the authors emphasize the discussion 
on the effects of temperature gradient at the wall and 
fluid temperature which may vary over time at the 
entrance as it is the case in systems cited previously. 
In order to obtain a more realistic behavior of engine 
with oscillating flow, as mentioned previously, the influ¬ 
ence of temporal variations of the inlet fluid temperature 
and its phase lag with pressure is particularly discussed 
in this study. 


2. Overview and model parameters 

2.1. One-dimensional problem 

The considered heat exchangers are currently made of 
many parallel channels so the basic geometry consists of 
two plane-parallel plates with an inter-plate space e. The 
channel height (see Fig. 1) is considered as infinite in the 
z direction and the length is limited to l 0 in the x direction. 


/~ V 
/ / 

/ / 

/ / 

/ / 


/~ 7 
/ / 

/ / 

/ / 

/ / 



Fig. 1. Basic studied geometry, oscillating flow in the x direction. 


Solid walls may have a temperature gradient along x which 
generates thermal coupling phenomena between the oscillat¬ 
ing fluid and the walls. 

The oscillatory motion considered takes place along the x 
direction and has an amplitude “D”. In Ref. [17] the authors 
show that general two-dimensional equations for the flow 
and temperatures can be reduced to a one-dimensional phys¬ 
ical model in the x direction as the system is supposed infin¬ 
ite in the z direction. The method consists in averaging 
fluid velocity and temperature in the y direction, 
( u, T g ) = (l/2c) («, T g )dy so the model [17] describes 

the averaged values of physical parameters (p , T, p. u, etc.) 
within any (y, z) section. Thereby, it is necessary to introduce 
friction factor and heat transfer coefficients for the wall 
boundary values. In the model, the thermo-physical proper¬ 
ties of fluid and material are also assumed to be temperature 
independent. All parameters introduced during the non- 
dimensional treatment are listed in Appendix 1. Spatio- 
temporal evolutions of fluid temperature depend both on 
fluid/wall thermal transfer coefficients and thermal condi¬ 
tions at plates’ extremities: end of solid wall, mixing of 
fluid, connected devices, and existing heat exchangers. 
All these parameters are functions of the relative particle 
displacement A L and of the plate length l 0 . As explained, 
we have to introduce in the ID model different correla¬ 
tions for channel geometry and to take into account pres¬ 
sure losses as well as oscillatory heat transfer boundary 
conditions at the wall. 

The wall friction factor Cf and the fluid to plate thermal 
transfer coefficient h are both expressed by means of a com¬ 
plex form in order to introduce two phases: the phase lag 
between fluid velocity and wall stress and the phase lag 
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between thermal flux and fluid/wall temperature difference 
that may occur in oscillating flow mode [19]: 


Cf = — 


du \ 

2/U a r| e / 2 
p \u\u 


and 


Nu = 



-d h 


ar 



( 2 ) 


(3) 


These formulations are detailed in the previous paper [17]; 
parietal and mean values of velocity or temperature are com¬ 
plex, as usual in linear thermoacoustics. The forms 
u = uofye^' and li = Uoge^'“ are adopted, where the two 
functions/j,,g are depending on the canal geometry. Real 
and imaginary parts of Cf and Nu were calculated in [18] 
for rectangular channels and represented versus the 
Womersley number. 

In the correlations for Cf and Nu, the Reynolds number 
can conventionally be based on the hydraulic diameter d h 
( d h = 2e for two infinite planes) and on the axial velocity 
magnitude: 


dhP gH 
Mg 


(4) 


Equation of the fluid energy (basis of the present study): 


VJI , . Nu, 1 = 


Nun 

4iNWiJW + AhC ° St dF + 


APrWl Tg 


1 

1 4 PrWL 


Nu,R\T t + Nuj\ 


dr 


(9) 


A usual simplification (see Section 2.3) consists in imposing 
a constant linear temperature at the wall: 



( 10 ) 


where T t and T_ i are imposed wall temperature at extremi¬ 
ties (see Fig. 2) 

Boundary conditions for the flow are taken as: 


x* = -l, U* = cos f, y;(r- 1) —,7 S , (11) 

X* = -1, Tf = cos f, r;(f,i) = r g+1 ( 12 ) 


Eqs. (11) and (12) postulate that when fluid enters the canal 
at x* = — 1 (left side in Fig. 2) or at x* = 1 (right side in 
Fig. 2) it is at a uniform temperature but this temperature 
may depend over time (Section 2.4) 

To solve the energy equation, the x Eulerian variable is 
replaced by the Lagrangian variable: 

? = x* — Asin f (13) 


Hence, in the non-dimensional form (stared notation), 
boundary conditions at the walls for velocities and thermal 
flux can be expressed by these more appropriate forms: 





(5) 



With the assumption of an incompressible fluid (p* = 1) 
in an oscillating flow mode between two plates, the 
characteristic equations of the problem are detailed in 
Ref. [17]: 

Fluid velocity is uniform along the canal and time 
dependent: 

2 u 

u — -= cos t (7) 

ooD v ' 

Equation for the pressure drop (not analyzed in this work): 


4 A l dp* 
~Ml dx* 



3[c f ]fo \ 
8 ^ m ) 


sin t * 


R[C r ] Re 

8WL 


cost* 


( 8 ) 


Then, the following two parameters are introduced (see 
other numbers’ definitions in Appendix 1): 

4 

Nu n +4PrW 2 om 
_ Nu r i 
^~Nun+4PrW 2 om 

Finally, Eq. (9) can be summarized by the following general 
form: 


(14) 

(15) 


g^ + /Tr; = £r;(? + Asint*) (16) 

The |3 parameter depends on the imaginary part of the 
heat transfer Nu, i, on the Prandtl and Womersley numbers 
and characterizes (Eq. (22)) the phase lag of the fluid to 
wall heat transfer. The A parameter is important when 
wall temperature is not uniform, so it characterizes the 
“regenerator” effect in that case, i.e. possibility for 
the wall to restitute to the fluid the heat deposited along 
the previous semi-period. 

The evolutions of the two parameters {X/A\f) and (3 as 
a functions of W om have been studied on Fig. 3 in 
Ref. [17]. The Prandtl number has the value Pr= 0.7 and 
the Nusselt numbers Nu Rl and Nu n are estimated by using 
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Fig. 2. Description of the oscillating flow in the x direction according to the ratio A\ . The scheme represents half an alternation of the fluid 
displacement. The central rectangular box is the volume between two plates of the heat exchanger. The fluid flows from the left (initial po¬ 
sition) to the right (final position). From top to bottom, Ar = 1 stands for heat exchanger where displacement of fluid is equal to exchanger 
length, Ar < 1 for long exchangers and /l i j > 1 for short exchangers. 


relations given in Appendix 1 [17]. For ILom varying from 
0 to 50 (A/Al) varies between 0.835 and 1 because the imag¬ 
inary part of Nusselt number has a week influence. /? quickly 
decreases when W om increases: since W om > 1.3 and /3 < 2, 
W om and |8 decrease asymptotically to 0. 

2.2. Case of a constant and uniform wall 
temperature T t = 1 

In most real cases, “entrance” fluid temperatures (and 
also the pressures), on the two extremities , T K \ for 
non-dimensional co-ordinates x* = — 1 and x* = 1, vary 
over time and are out of phase with respect to fluid displace¬ 
ment. In a first step, in order to facilitate understanding, we 
consider their values as constant and not time dependent; 
variations over time will be discussed lately. The arbitrary 
choice of a hot fluid temperature at the left inlet 
(x* — — 1) and a cold one at the right inlet (x* = 1) is 
made; the wall temperature is kept constant and uniform. 

For constant boundary conditions, the solution for gas 
temperature is quite simple: 

Ar rad (f*,?) = (T;(C?) - l)e-«'-a (17) 

where f* n corresponds to the fluid particle entrance time at 
x* = —1 and A7’ ra d(t*,f) = T*(f*,f) —T* is the fluid to 
wall temperature difference in Lagrangian formalism f at 
time t* and along the system. 


Thereby, Eq. (17) represents the history of the tempera¬ 
ture difference between the gas and the wall. If the consid¬ 
ered gas particle comes from the initial position xj, at time 
tg, and enters into the plates at time t* n , its position at any 
time t* is then given by: 

x* = xj + A l (1 + sinf*) (18) 


Eq. (17) is only valid between t* n and f* ut , the latter instants 
defining the times when the fluid particle remains within the 
system. f* n and f* ut are given by: 


t. = arcstn 


0 +*0) ^ 

) 


? out — 2 ? o + 2 tc t m 


(19) 

( 20 ) 


The temperature difference provided by Eq. (17) appears in¬ 
directly as a function of the fluid displacement along the 
plates by the value of the fluid where the non-dimensional 
particle position x* is masked by the Lagrangian variable f 
(see Eq. (13)). Moreover, due to thermal transfers between 
plates, an exponential decrease (/3 parameter) of the initial 
(in x* = — 1) wall to fluid thermal difference AT rad is 
observed, suggesting an optimal exchanger length may exist 
for a maximal heat transfer. 

In fact, the fluid particles located between the plates and 
initially located at xj) > — 1 (if time origin is tj = — (tt/2), 
and then u* = 0) never reach the position x* = — 1 because 
their displacement is on the right. On the other hand, 
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according to the particle displacement value A l, a part of 
these particles moves up to x* = +1 (see Fig. 2). 

To simplify and to clarify the discussion, we first con¬ 
sider the fluid particle that enters a well mixed and constant 
temperature tank on both sides of the canal. At t* = tg + n, 
the direction of the fluid particle changes so the movement 
occurs from x* = 1 to x* — — 1. In the model presented 
here, the fluid particles flowing back (cold fluid) into the 
exchanger at x* = 1 are considered to be at a temperature 
inferior of Ar rad(+ i) = —0.75 to that of wall, and those 
(hot fluid) flowing into the exchanger at x* = — 1 to be at 
a temperature AT rad( _ 1) = +0.75. In the case of the relative 
displacement of fluid particle which is less than unity 
(A l < 1), a part of fluid remains inside the exchanger forever 
and never flows out. Therefore, its temperature reaches the 
wall temperature exponentially over time (Eq. (17)). This 
fluid part forms a “fluid stopper” zone that does not contrib¬ 
ute to any global thermal exchange. 

Fig. 3 shows the modification of the fluid temperature 
profile along the plates for three exchanger geometries: short 
A l = 1 .5, with a length equal to fluid displacement A L = 1 
and long H L = 0.5. The results are presented considering 
the thermal parameter (3 = 0.5. Hot fluid enters the 


exchanger at constant temperature at x* = — 1 and cold fluid 
in the reversed flow at jc* = I. Straight lines represent, time 
after time, the fluid temperature profiles during the move¬ 
ment from left to right, and dotted lines, the profiles corre¬ 
sponding to the movement in the opposite direction. One 
can clearly observe the exponential decrease of the thermal 
discontinuity front occurring time after time during the fluid 
progression. In Fig. 5 of [17], it is shown that the smaller (3 is 
(for large W om values), the bigger the observed temperature 
front (since fluid entrance) is, i.e. gas/wall temperature dif¬ 
ference remains high hence indicating that thermal transfers 
are quite bad probably. This clearly indicates that the rate of 
the thermal front decreases, thus exchanger efficiency de¬ 
pends upon Hl- In fact, only fluid particles reaching both 
the exchanger extremities at x* = — 1 and x* = 1 are able 
to exchange heat with the wall. 

Here on Fig. 3, for a similar absolute fluid displacement, 
one notes that a too long exchanger (Hl = 0.5) has an inactive 
zone for both hot and cold fluid, respectively, between 
.r* = — 1 and 0 and x * = 0 and +1. On the contrary a too short 
exchanger ( A L = 1.5) does not permit a maximum heat ex¬ 
traction. In fact, “exit” temperatures atjr* = +1 and ’"return” 
temperatures at x* = — 1 are higher in such an exchanger. 




A, =1.5 



x* 

Fig. 3. Histories of the non-dimensional temperature profile along the exchanger [—1, +1] and of the relative fluid displacement, case of a con¬ 
stant wall temperature A7| orlg — 0 (horizontal thick line). Hot fluid inlet AT rad( _i) = 0.75 on the left on the first half alternation, and cold fluid 
inlet on the right AT rad( i) = —0.75 on the reciprocating flow. Full lines represent a flow period from + = —1 to x* = +1, dotted lines in the 

other way — /3 = 0.5 and Pr = 0.7. A\ — I stands for heat exchanger where displacement of fluid is equal to exchanger length, A\ . I for long 

exchangers and > 1 for short exchangers. 
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After some tedious calculus [17], the global heat quantity , t t A 

exchanged between fluid and plates during half period (left Ar rad [t ,x ,x 0 ) + ^\Ti ong cos(t <p) 
to right for example) can be expressed by: 

Grad = PC t S f u.T 0 *T a H- 1) -JU? l tf) (21) 

1 With the phase tp depending on (3: 


Ar rad (f\ — 1) + -AT long c 


»(£-*))' 




(24) 


where F t is defined by: 


Fiifi) 


1 + 8/? 2 - e- 2(37c 
1+4/3 2 


( 22 ) 


Since the maximum heat quantity that could be exchanged 
by the gas along its displacement through the exchanger is 
2max i — pCgDS , p ass r 0 Ar rad( _i ) , the function (l/2)Fi 
represents the exchanger effectiveness value: 
E = (2rad/2max/) = (l/2)Fi. Fig. 4 shows this effective¬ 
ness as a function of (3 parameter: the maximal value of 
100% is reached for f3 > 2 (Pr = 0.7). When W om tends to in¬ 
finity, /? and the heat quantity exchanged tend both to 0: 
thereby at high frequencies, thermal transfers may vanish. 
That means the inter-plates distance «e» must be reduced 
to ensure a low value of W om and an adequate thermal ex¬ 
change (since the thermal boundary layer diminishes with 
frequency, then radial exchanges become less important 
than axial diffusivity). 


2.3. Case of a linear and time constant wall 
temperature profile 


The wall temperature profile is now given by a linear 
equation: 


r p* _ r p* 

(*•)= ■ _ -v 


+ 


Tt + Tf, 


(23) 


Thermal transfers phenomena are more complicated when 
a temperature gradient along the plates’ walls is imposed: 
in the case of particles flowing into the exchanger and passing 
on — x* = —1, the expression of the fluid to wall local tem¬ 
perature difference for x* and t* becomes ([17], Eq. (52)): 


cos (p = 


v / i+^ 


sin (p = 


Vi + P 2 


(25) 


And the time interval when a fluid particle remains in the 
canal is given by: 

(l +*o) 


t- = arcsin 


1 


Al 


; and C = - 4 


(26) 


Let us now express the wall extremities temperature differ¬ 
ence (constant) by: 

AT long = T* - r j (27) 

According to the wall temperature profile (Eq. (23)), the 
fluid to wall temperature difference in Lagragian formula¬ 
tion is: 


ATU4?) = r;(t*,?)-r,V) (28) 

It is always depending on the position f = x* — Asin t* of the 
fluid particle at the time t*, but now local wall temperature 
varies along a*. 

Eq. (24) is valid only for particles never exiting at x* = 1 
and over times t*e [t* n , t* u J. So the two temperature differ¬ 
ences AT rad( _n, AT rad (i) and the wall extremities tempera¬ 
ture difference A7] ong both contribute to the variation of 
the fluid temperature profile; moreover an exponential 
decrease of parameter 0 still occurs. 

For fluid particles flowing back into the system and pass¬ 
ing through the inlet section at x* = 1 (\T ra $(t*, 1) = —0.25 
on Fig. 5), the following relation must be applied: 

Ar rad (r*,x* : x*) = —^A7) ong ^cos(t* n / - ? ))e“' i( ''“4> 

- cos (f - <p)) + AT rad (f*, l)e-«'*-4> (29) 



Fig. 4. Exchanger effectiveness representation (corresponding to 
F i/2) as function of parameter (3. 


With the time interval: 


- 1 


( 1-4 


A, 


and tl, = rr — f 


(30) 


Fig. 5 shows temperature profiles histories obtained 
along the plates within the fluid for the same three relative 
displacements Al and assuming the fluid inlet temperature 
is the same of Fig. 3 (so AT rad (?*, 1) = —0.25 and 
AT rad (t*, —1) = +0.25 in order to respect this hypothesis). 
Whatever Al and (3 values, the existence of a thermal gradi¬ 
ent along the walls induces the “regenerator” effect for ther¬ 
mal exchanges as explained in the first part of this work. A 
part of the heat quantity exchanged between gas and wall 
may be restituted in the second half period; for example, 
this situation must be avoided in the case of a cooler 
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Fig. 5. Histories of the non-dimensional temperature profile along the exchanger [—1, +1] and of the relative fluid displacement, case of a 
linear profile of the wall temperature A7\ ong = 0.5 (thick line). Hot fluid inlet A7 nK |, |> — +0.25 on the left on the first half alternation, and 
cold fluid inlet on the right Ar rad( i, = —0.25 on the reciprocating flow. Full lines represent a flow period from** = — 1 to x* = + 1. dotted lines 
in the other way — ft = 0.5 and Pr = 0.7. Al 1 stands for heat exchanger where displacement of fluid is equal to exchanger length, Al < I 
for long exchangers and Al > 1 for short exchangers. 


exchanger. After the flow direction change, the fluid parti¬ 
cles that are colder than the wall (located under the thick 
line of wall temperature in Fig. 5) have a reversed heat ex¬ 
change because they move along a positive thermal gradient 
hence the fluid is heated back. That is the reason why this 
process can be qualified as a wall “regenerator effect”. At 
the same time, particles that are still hotter than the wall 
cool anymore. When the thermal gradient through the wall 
is negative, the regenerator effect is added to the exponential 
decrease of the “entrance” temperature gradient. In the par¬ 
ticular case where AT ra(1 /, n is strictly symmetric to Ar ra d(_i) 
with respect to wall temperature (case of Fig. 5), the global 
thermal balance becomes null and the exchanger operates as 
a regenerator actually; on the contrary if Ar rad ( +1 ) and 
AT ra( ](_r) are of the same sign, thermal exchanges are 
doubled. 

Calculus of the exchanged heat quantities leads to distin¬ 
guish two cases depending upon the exchanger length: short 
(Al >1) and adapted or long (Al < 1): 

• case of adapted or long exchangers (Al < 1), with 
AT rad ( +1 ) = 0 

Qa l < 1 = Grad = PCgSpass^r 0 Ar rad (_i)+l L F 1 (/?) (31) 


• case of short exchangers (Al >1), with AT rad ( + i) = 0 

G+>i = P c gSpass^ 7 ’oA 7 ’ rad( _ 1 )(Gii +G21) 

— pc g S pass |7'o-Ari ong (Gi2 — Q P i) (32) 

Note: if AT rad / +1 ) =£0, one must add a supplementary 
term similarly to Ar rad (_i) 

Terms Git> Gi2> G2n-->G P i are demonstrated and calcu¬ 
lated in Ref. [17] (Eqs. (62), (63) and (65)). They correspond 
to integral heat quantities exchanged for particular volume 
of fluid during each half period (fluid on the left or on the 
right or inside of the plate system at initial time, Fig. 2). 

Once again the evolutions of the right member terms of 
Eq. (32) still have been represented on Figs. 8 and 9 of 
Ref. [17], The term (G11 + G21) represents a degradation fac¬ 
tor of the efficiency (with respect to the case of an exchanger 
of sufficient length) because the ratio (G11 + Qii)/AlF\ is 
always inferior to one. Its value is lower especially when 
Al increases and (3 decreases (large W om ): therefore, in the 
case of a too short exchanger, the thermal exchange is not 
efficient. Otherwise the last term (G12 — Qpi) remains 
always negative and decreases when Al increases. The re¬ 
generator effect and the exchanged-restituted energy parts 
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are characterized by this second term in Eq. (32). The global 
sign of the second term of Eq. (32) depends on the sign of the 
temperature difference Ari ong . 


2.4. Fluid inlet temperature depending on time 

The functioning of machines like Stirling, Vuilleumier, 
Gifford-Mc Mahon, Pulse Tube Refrigerator or thermo- 
acoustical systems is often quite different from conditions 
previously described; it may be characterized by: 

- the use of a compressible fluid, 

- the flow is usually turbulent, 

- fluid and materials physical properties may strongly 
vary (at cryogenic temperatures for example), 

- periodic fluctuations and phase lag fluctuations of 
thermodynamic quantities: pressure, temperature, dis¬ 
placement, velocity, etc. 

- the amplitude of these fluctuations that may not match 
the linear thermoacoustic theory. 

When Ar rad( _i) and Ar rad(+ i) vary over time sinusoidally 
(like in Stirling machines for example), the previous expres¬ 
sion of F i (/?) appears as a function of the phase lag between 
the temperature and the inlet velocity. Therefore both its 
integration and that of initial Eq. (16) become more com¬ 
plicated actually. 

Let now pay attention onto the modifications that must be 
applied to previous results in the specific following case: the 
fluid inlet temperature at the device extremities varies over 
time and its phase lag with the fluid velocity variations is 
tpj. This temperature variation is supposed of the following 
periodic form: 



Note: velocity is still u* = cos t* and in Lagrangian 
formulation r ; * n and x* depend both on the initial particle 
position Xq. 

Fig. 6 illustrates the thermal front of the fluid, along the 
canal at different instants over a period and for a particular 
value of the temperature phase tp T = (it/2). For simplic¬ 
ity, the case of a null ATi ong is chosen; otherwise the 
regenerator effect still exists and is not influenced by 
the variation of the fluid inlet temperature. One can note 
the important difference existing in comparison with the 
previous graphs of Fig. 3. The graphic representation is 
the same: full lines stand for each instant when fluid flows 
from x* = — 1 to x* = +1, and dotted lines in the other 
way. Horizontal arrows symbolize the fluid displacement; 
vertical arrows represent the variation of the fluid inlet 
temperature amplitude at the same time. One clearly ob¬ 
serves the combination at side x* = — 1 of the diminution 
of the inlet temperature during the displacement of the 
fluid toward the left and the exponential diminution of 



Fig. 6. Temperature profiles histories when the inlet fluid tempera¬ 
ture varies over time for a heat exchanger where displacement of 
fluid is equal to exchanger length Al = 1. Case of a constant 
wall temperature A7i ong = 0 (horizontal thick line). Hot fluid inlet 
amplitude Ar rad( _!, — +0.25 and phase <p 7 , = ~/2 on the left for 
the first half alternation, and cold fluid inlet on the right 
Ar radll) = -0.25 and phase <p T ,, — ~/2 for the reciprocating 
flow. Full lines represent a flow period from x* = — 1 to x* = +1, 
dotted lines in the other way. Vertical arrows indicate the time inlet 
temperature evolution of the fluid —0 = 0.5 and Pr = 0.7. 

the hot thermal front due to thermal transfer with the 
wall. After the alternation of the fluid, the temperature 
of the cold fluid drops as long as the fluid is entering at 
x* = +1 while the heat transfer between the wall and 
the hot fluid is going on. When evaluating the heat quan¬ 
tities, only the terms related to the previously constant 
temperature difference Ar rad( _ 1) have to be substituted 
by Eq. (33), thermal exchange phenomenon due to the 
gradient ATi ong being not affected. 

Fig. 7 shows the new term related to exchanger effi¬ 
ciency that replaces the value A^F j(/3) in Eq. (31) for 
long or adapted exchangers; it is reported in a 3D graphic 
as function of the phase lag cp T and phase thermal transfer 
parameter 0: A section of the graph in the plane 



Fig. 7. Efficiency factor AlF\ as a function of 0 and of the phase 
lag tp T between the gas temperature and the fluid velocity. 
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Fig. 8. Corrective term of heat quantities (Q n + Cbi ) /{A\J'\ (3)) 
as a function of Al (Al > 1 and <pr = 7t/4). 


<p T = const shows that the evolution is similar to that of 
efficiency in Fig. 4, however, according to <p x the maximal 
value is very different, and even null if <p T = 0. still 
has to be chosen as high as possible, in order to limit the 
Womersley number. 


A h F\ = A h 


cos 2 <p T e~ 2 ^ n [Cl — C2] 


C 3 


(34) 


With: 


Cl = (1 +4/3 2 ) (^1 + e 2 ^ IC (2 l S 3 Tu + 2/?Tc— l)jcos(p T 

C2 = /3(4(1 + /3 2 ) ((l + 8/3 2 )e 2#iTC — l) 

+ (l + 4/3 2 ) (e 2 ^ — l)sin <p T ) 

C3 = 8/3 2 (l + 5/3 2 + 4|S 4 ) (sin 2 ^ - l) 

The correction term Qn + Q 2 \ °f heat exchange efficiency 
for short exchangers (Eq. (32)) has also to be modified. 
Only numerical values can be obtained: an example of the 


relative value compared to AlF\ ((3) for a temperature phase 
lag of cp T = 45° is given in Fig. 8. In comparison with Fig. 8 
of [17], the values are quite lower in comparison with the 
case of a long exchanger of constant parietal temperature. 


3. Average density of enthalpy flux 


A good criterion for evaluating thermal exchanges be¬ 
tween fluid and wall in an exchanger operating in a periodic 
flow regime consists in observing the variations of the local 
and time averaged enthalpy flux h = p g c g uT and 
H = (w/2tc) j' 1 ™ h(t)dt. In dimensionless form and for 
the first order, i.e. with a constant density, h* and its average 
value H* can be expressed as following: 


, , , 2 h 

h * = p*«*r = - - 

s toDp 0 T 0 c„ 


and H* = p*(u*T\ 


2 H 

toDp 0 T o c g 


(35) 


The goal with a heat exchanger is to maximize the variations 
of the average enthalpy flux between the two device extrem¬ 
ities (an increase in a heater, a decrease in a cooler). The ID 
averaged model developed here makes easy the study of heat 
exchangers in oscillating flow mode, but conditions of real 
functioning being complicated, the understanding of phe¬ 
nomena remains sometimes fragmentary: boundary condi¬ 
tions, phase lags, temperature profiles, etc. 

Fig. 9 shows axial variations of the average enthalpy H* 
for different kinds of exchangers operating with or without 
a thermal gradient at the wall (same values to those of 
Fig. 5). Variations of H* are more significant if there is an in¬ 
crease of j8, taken to 0.5 for example. For an exchanger in 
which length represents twice the displacement value, half 
the device does not contribute to heat exchanges: indeed, in 
such a case, the global thermal behavior is similar to that 
found for Al = 1. With a too short exchanger, global 


p = 0.5 


P = 0.5 < P n -i) = rr / 2 <Pn+i) = 7 Z, 




Fig. 9. Axial variation of the average enthalpy H* — p* (u*T*) for different kind A l of exchangers of uniform wall temperature (/? = 0.5), full lines 
AFiong = 0, dotted lines Ari ong = 0.5 left: constant inlet temperatures, right: variable over time inlet temperatures with phase lag 90° on both sides. 
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variations of H* are not very different, however, temperature 
levels need to be higher and the axial profile form is very dif¬ 
ferent. Eventually, the thermal gradient in the wall (negative 
in the above example) induces a decrease of H* in all cases. 
Influence of phases of inlet temperatures is more difficult to 
explain: in the chosen example, it is clear that enthalpy level 
diminishes because average inlet temperatures also diminish. 

4. Conclusions 

It has been demonstrated for oscillating fluid flows that 
classical formulations of thermal heat transfer evaluations, 
and especially those corresponding to steady fluid flows re¬ 
gime, cannot be employed. This is a consequence of the 
phase lag existing between all the physical variables (pres¬ 
sures, displacements, velocities, temperatures, etc.) and the 
energy transfers. This is also a consequence of the particular 
interaction between the wall profile temperature and the 
oscillation of the fluid. 

Two important parameters (in addition to Prandtl num¬ 
ber) characterize oscillating flow thermal transfers: Al, the 
ratio between fluid displacement and heat exchanger length. 


and the Womersley number W om depending on the oscillat¬ 
ing frequency and the hydraulic radius. The model studied in 
this paper allows the calculation of the thermal effectiveness 
(F\i 2 ) of heat exchanger with incompressible oscillating 
flows; this is a first step because compressible flows are 
too complicated for our current knowledge. Corrections 
functions are given for the case of short heat exchangers 
Al > 1, when a longitudinal gradient exists ATiong^O, or 
when fluid inlet temperatures are time dependent. 

Meanwhile, conclusions are still interesting: 

- it is unprofitable to use too large or too short heat ex¬ 
changers, better performances are obtained when there 
is adaptation between length and fluid displacement; 

- when a temperature gradient is present along the wall 
of the exchanger, the efficiency can be dramatically 
diminished by a “regenerator effect". 

The last remark is particularly important for designers of 
heat exchangers, Womersley numbers superior to 1.3, corre¬ 
sponding to /3 < 2 (case Pr = 0.7), are prohibited because the 
thermal efficiency dramatically vanishes, so hydraulic radius 
have to be adjusted according to the operating frequency. 


Appendix 1. Complete expressions of exchange and friction coefficient 

Plates: 

Nu _ 4jPrW^nh(W om VJP?) 8jW 0 >nh(W omV //) 

W om \fjPr — tanh(IT om \fjPt ')’ ‘ W omS fj - tanh(W omV {/) 

Tubes: 


Nu = —2W„ 


/2VFr MW^Vp-rf'i) 

Jo{W cm VPtf/ 2 ) - 1 


= 4W om / /2 


W om fF 


uw om f/ 2 )y' 

A (WW 3/2 V 


Table 1 

Main dimensionless numbers characterizing oscillating fluid flows (from ID model) 


Numbers definitions Name 


Q 1- 
II 

J 

Non-dimensional ratio 

'0 


w /Twf 

W -=\J M 

Womersley/Stokes 

(i)D 

M m - _ 

yPo/po 

Mach 

PrM 

Prandtl 

CfR e 

Darcy Product 

aj 2 eft[h\ 

NUri — 

& 

Real Nusselt 

Nu n = 2em 

Imaginary Nusselt 


Signification 

Ratio of total acoustical displacement to axial length 
Important for extremities effects if D is large in front of /q 


Ratio of hydraulic radius to viscous diffusion thickness 


Ratio of acoustical velocity to sound velocity 

Acoustic effects, variable compressibility effects to be considered 

If M m >0.1 non-linearity of the phenomenon 

Ratio of viscous and thermal diffusion length 

Characterizes local frictions on walls 

Characterizes the external thermal transfer 


Characterizes the internal thermal transfer into walls 
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